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In the framework of the paraxial and of the slowly varying envelope approximations, with reference 
to a normally dispersive medium or to vacuum, the electromagnetic field is given as a continuous 
quantum superposition of non-dispersive and non-diffracting wave-packets (namely X- waves). En- 
tangled states as pairs of elementary excitations traveling at (approximately) the same velocity are 
found in optical parametric amplification. 

o 

O . INTRODUCTION 

<N ■ 

The quantum description of diffracting non-monochromatic beams has been the subject of intense research (see 
fc^H ' e.g. references 

QH 111113) Various approaches have been proposed in order to provide a consistent quantum 
picture of nonlinear optical processes and related phenomena, like entanglement and photon localization. Several 
' important issues have been discussed in the literature, as equal-time or equal-space commutation rules, or paraxial 
approximations, and different "schemes" have been formulated. [U 0, E3, El EH Most of them are based on the 
second quantization of the relevant classical equations, or on a three dimensional plane wave expansion, however their 
employment in specific problems is far from being trivial. In this Letter I describe a new formulation for 3D optical 
| field quantization, and discuss its application to optical parametric amplification. 

All the machinery is based on the paraxial and slowly varying envelope approximation (SVEA), with reference to 
vacuum and to a normally dispersive medium. However a general theory can be envisaged. The basic element is 
recognizing that a traveling pulsed beam can be expressed as a continuous set of harmonic oscillators, each weighted 
by a rigidly moving 3D wave-packet. It corresponds to a non-dispersive non-diffracting solution of the relevant 
classical wave equation, with velocity v. This is the so-called X-wave transform, or X-wave expansion, as outlined 
below. [Tllll2| Then the quantization comes naturally, following standard approach, using continuous-mode operators. 

US 03 

X-waves are 3D invariant pulsed-beams trav eling wi thout diffraction and dispersion, also known to as "progressive 
undistorted waves" (see e.g. the recent reviews lldlld and references therein). Their appearance in nonlinear optical 
processes has been experimentally investigated in |17| . 

For free propagation, it turns out that the quasi-monochromatic beam, with carrier angular frequency oj, behaves 
like a quantum gas of particles which, in vacuo, have mass m satisfying the relation mc 2 = Tllo. The rigidly moving 
. "eigenmodes" (the X-waves), while not having a finite norm (like plane waves), do have a strong degree of localization, 
thus addressing the issue of spatio-temporal photon confinement (see e.g. and references therein). In a simple 
way, this approach provides the quantum description of the 3D evolution of a classical, paraxial and slowly varying, 
beam. Furthermore classical infinite-energy X-waves, here described as coherent states, do have finite expectation 
values for energy in the quantized picture. 

DERIVATION OF THE PARAXIAL EQUATIONS 

In this article a new quantum-mechanical description of pulsed light beams is introduced. In order to avoid as much 
as possible unnecessary complexities and involved formalism, paraxial slowly varying beams are considered. Given 
the fact that progressive undistorted waves are well-studied solutions of the wave equation, a much more general 
theory can be envisaged. Typical experiments in nonlinear optics arc described in the framework of the considered 
approximations and, for similar reasons, the treatment is limited to circularly symmetric beams. 

The basic model can be obtained by the scalar 3D wave equation for the electromagnetic field E in a medium with 
refractive index n: 



(1) 
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By letting E — Re[Aexp(—iujt + ikz)]/ yje^n 2 /2, with k = ujn/c, in the paraxial approximation the usual Foch- 
Leontovich equation can be derived: 

,8A n w9i _ 9 . 1 d 2 A n 
2 ^ +2 ^ +V ^-?^ = °- (2) 

Equations like J5J) are typically adopted in nonlinear optics, and interpreted as a propagation problem along z. If 
A, as a function of t, is slowly varying, they can be equivalently cast as a time-evolution problem; the latter is the 
formulation adopted here in order to be as much as possible consistent with standard quantum mechanics. With this 
aim, introducing £ = z — ct and r = t, Eq. is rewritten as 

If the SVEA is applicable, the right-hand side of the previous equation is negligible and, coming back to the original 
variables (z,t), Eq. can be written as 

8A ,8A uj" OA u)' o 

l ^ + ^-T^ + 2fc V ^ = ' (4) 

while being u>' = c and u>" — c 2 /oj. Eq can be treated as an evolution problem with respect to time t. 

The very same procedure can be repeated in a dispersive medium, as done in Ref. 0, while retaining terms up to the 
second order dispersion. The model still holds with uj' and uj" the derivatives of the dispersion relation uj = ui(k). 
Hence, since it can describe vacuum as well as a dispersive material, [3(| and it is typically used with reference to 
experiments with pulsed-diffracting beams, it will be the starting point for the quantization procedure. Without loss 
of generality the envelope field A has been normalized so that the time averaged electromagnetic energy in the whole 
3D space is 

8 = J jj \A\ 2 dxdydz, (5) 
i.e. any type of box-normalization is not necessary. 



THE X-WAVE TRANSFORM AND X-WAVE EXPANSION 



Assuming radial symmetry hereafter (r = y ' x 2 + y 2 ) , the general solution for the initial value problem is expressed 
by the Fourier-Bessel integral 



OO pOQ 

ik z Z—iClt , 



A= / fc_LJo(fc±r)5(fc ± ,fc z )e^ z - lS "dfc ± dfc z , (6) 

J -oo JO 

with Z = z — ui't. Eq. © furnishes the field at instant t, given its spectrum S at t = 0, with transversal and 
longitudinal wave-numbers k± and fc z , and ft — —uj"k 2 /2 + uj'k 2 _/2k. With the following change of variables (a > 0) 



k± = otyJuj"k/u)' 

I have, after some manipulations, 



k z =a-v/u>", (?) 



OO 

A= I e-i&r^vfaZ-vQda, (8) 

while being 



DO 



II 



9 v (r,Z-vt)= I X{a,v)J {\J l ^ar)e l{a -^ ){z - vt) da. (9) 



Hence A can be expressed as a superposition of self-invariant beams (i.e. waves traveling along the z— direction without 
dispersion and diffraction) with different velocities. The quantity X(a, v) — ^rS(aJ ^^,ol — -^t) corresponds to the 
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so-called X-wave transform. [IlJ X- waves are defined as the propagation invariant solutions of the wave equation J3J, 
for more details see Ref. and references therein. 

X(a, v) can be expanded by Laguerre functions (2x)e~ x with respect to a: 

X(a,v) = -£%L C p (v)f p (a) 

(10) 

/» = V / ^4+T)( aA )4 1) (2«A)e-« A , 

where A is a reference length that is related to the spatial extension of the beam (the specific form of f p is chosen for 
later convenience). A is written as a superposition of basis orthogonal X-waves (given by 10 with X(a,v) = f p (a)): 



A = X™ =0 CpMe-^VpfaZ -vt)dv. (11) 

J —OO 

The following relation holds by construction 

< ^(r, Z - ut)\ip;(r, Z - vt) >= 6 pq 6(u - v), (12) 

where < f\g > denotes the integral over x, y, z of /* g. Eq. (|12|) reveals that the X-waves have infinite norm, like the 
plane waves typically adopted for field quantization. 



QUANTIZATION 

The classical energy of the pulsed beam, Eq. (J5J , is given by 

/oc 
\C P (v)\ 2 dv, (13) 

while the field in can be rewritten as 

/oo 
C p (v,t)r p (r,Z-vt)dv, (14) 
-oo 

being 

^{v,t) = -iu p {v)C p {v,t) (15) 

with Up — v 2 /2lo" . Hence A is the integral sum of harmonic oscillators, with complex amplitudes C p , each associated 
to a traveling invariant beam (corresponding to the usual cavity mode, with the difference that it is rigidly moving) . 
Introducing the real- valued "position" Q p (v) and "momentum" P p (v), by letting 

C p (v,t) = -^=[u Jp (v)Q p (v,t)+iP p (v,t)}, (16) 



£ is written as the sum of the time-independent classical energies of each oscillator 
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f = S^ / \[P p {v,t) 2 +io p (v) 2 Q p (v,t) 2 }dv, (17) 



uj p (v) is the resonant frequency of each "mode". 

Following standard approach to quantization, |2(j l2ll | it is natural to write the Hamiltonian like (omitting extrema 
for the integral and the S symbol) 



H = J %uj p {v)a p {v)a p {v)dv = E p J ^-a p (v)a p (v)dv, (18) 

with \a q (u), a p (v)] — 5 pq S(u — v). Note that the zero-point energy has been omitted in (|18f) . as a standard renormal- 
ization in quantum field theory. |2d| 
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Some comments are needed at this point. First of all I observe that l|18|) corresponds to a set of freely moving 
particles, each contributing to the energy by the kinetic term, with mass defined by 

m=A (19) 

Thus the freely propagating pulsed beam is represented like a quantum fluid, or gas.j^ The eigenstates are given by 
\l,n,v >= [a|(u)] n |0 >, denoting n particles with velocity v in the traveling "mode" (or fundamental X-wave) ipf, 
and energy nrav 2 /2. Then, it is to note that, differently from standard heuristic approaches to field quantization, 
the Hamiltonian directly corresponds to the energy of the beam, and not to the second quantized form, i.e. to the 
Hamiltonian functional giving Eq. Q. 

The optical field amplitude operator is given by 



A — Sp 



: ipp [z - (u/ + v)t, r] huj p (v)a p {v)dv. (20) 



The mean value of A and of the energy density, in the state with one elementary excitation on the /—fundamental 
X-wave at velocity v, are: < 1, 1, v\A\l, 1, v >= and < I, l,v\A^A\l, l,v >= Tiuji(v)\i{j^ (Z - vt,r)\ 2 . 
Coherent states can we constructed as well. For example, taking (with obvious notation) 



|<>;(")| 2 ai(v) n 

|a,(«)>=e— ^-E~ ^f|I,n,i;>, (21) 



with ai(v) an arbitrary complex number, I have a state that has an expectation value for the optical field profile given 
by the fundamental X-wave ipi(Z — vt, r) and corresponds to a classical X-wave. |3l| However there is a remarkable 
difference: the expectation value of the energy is huii(v)\ai{v)\ 2 . While this state has infinite energy in a classical 
context, in a quantum mechanical picture X-waves do have finite energy. 

On the other hand it is well known, from quantum field theory, that \l,n,v > is not normalizable, and a, \(v) 
must be intended as operator valued distributions. 20] This is due to the fact that a state with a definite v is an 
idealization, like an elementary particle with a definite momentum k. Any experimental apparatus implies a spread, 
and an indeterminacy in the velocity (corresponding to a normalized state obtained as a continuous superposition of 
\l,n,v >). 

By particularizing Eq. (|19|l to the case of vacuum, using the previously mentioned relations, it is found 

mc 2 = Tiu), (22) 

which shows that the effective mass of the quantum gas particles is determined by the carrier frequency of the beam 
via the Einstein relation. 

It is also notable that, in vacuo (a/ = c), superluminal components are encompassed in Eq. (|20|l . But this is not 
surprising. Indeed the superluminal propagation of X-waves does not lead to paradoxes, since they are stationary 
solutions (i.e. they fill all the space, like plane waves) of the wave equations (see also 15, 23]). This can be shown to 
be perfectly consistent with the principles of standard special relativity, in the framework of the covariant Maxwell 
equations. 24] Furthermore it is well known that superluminal propagation has a not zero probability, and hence it is 
predicted, in standard quan tum electrodynamics. |25j Single-photon superluminal effects have been experimentally 
investigated in Ref. 1261 [32| Before considering an application of this formalism I want to underline that a completely 
different wavelet quantization of relativistic fields has been previously developed by Kaiser. |27j 



OPTICAL PARAMETRIC AMPLIFICATION 

I consider the phase-matched co-directional optical parametric amplification (OPA) of two frequencies, u>x and 
u>2, attainable in quadratic, or cubic, nonlinear media. For simplicity the pump beam is treated as a constant term in 
the Hamiltonian. A quantized field is associated to lo\ and u>2, with parameters oj'-^ 2 , ki t 2 — ^1,2^1.2/2. lo" is taken 
equal for both of them, and io' x 7^ u>' 2 . A\^ denote the corresponding optical quantum fields: 

Ax = Ep/e'^Vp.it 2 ~ ( w i +v)t,r] y / huj p (v)a p (v)dv ^ 
A 2 = ^ P J e~~ l ~t7? ip% >2 [z - (uj' 2 + v)t,r] y / huj p (v)b p (v)dv. 
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The basis X- waves are denoted by ip"j in order to indicate that they have to be calculated after Eq. (|10f> with 
parameters (co'j and kj) of ujj; the corresponding spectra are denoted by f P j(a) (j — 1,2). 

The total Hamiltonian is given by the sum of the free terms and the interaction Hamiltonian Hj: 

H = S p J huj p (v)(ala p + b p b p )dv + Hi. (24) 

Hj can be determined by the classical counterpart given by Ej — \ < Ai\A 2 > +c.c. (c.c. denotes complex conjuga- 
tion). After some (lengthy but straightforward) manipulations (h.c. is hermitian conjugate): 



Hj = hZ pq / / dudv Xpq (u + «)y / ^>H(^e^( u ' I ') t ap ( u )6t( v ) + hx . (25 ) 



where F(u, v) = (u 2 + v 2 )/2uj" + (u — v + — u)' 2 )(v — pu)/uj", with p = (k\uj' 2 jk 2 Lo']) x ^ 2 . The interaction term Xpq{ u ) 
is 



= ^X^A^^^jr^W, (26) 
v y k\ki (l + pja/' (l + p)a/' 

with 9{v) the unit step function. 

By applying standard time-dependent perturbation theory, assuming as the initial state |0 > at t = 0, at first order 
it is 

\^(t) >= H pq f J e iK ( u ^% q (u + v)G(u, v, t)y/Lj p (v)u q {u)al(u)bl(v)dudv\0 > , (27) 

with 

(u — v+uj'^ —uj' 2 ){v— pu) 



j\ ya, vj - - 2{l+p)w" 

„, f\ — 2sm[g(u,v)t/2] 



g(u,v) = + -^- pu) (28) 



The state given in Eq. (|27|l is a continuous-variables entangled superposition of particles at the two generated 
frequencies, traveling with different velocities. Indeed the weight function is not separable with respect to u and v 
(see e.g. Ref. flil and references therein). 

The transition probability to the two particles state a p (u)b q (v)\Q > is given by 

^ , % , . / \ i , m9 sin\q(u, v)t/2] 2 9 . „ N 

V pq (t,u 7 v) =u; p (v)cj q (u)\x Pq (u + v)\ 2 (29) 

As t — > oo, V pq tends to a Dirac-delta peaked at g(u,v) = 0. Entangled particles, with velocities u and y , are 
associated, for large propagation distances, to a point in the plane (u,v), lying on the parabola g(u,v) — 0. |3.1 | Of 
particular relevance is the low velocity region. 

Proceeding as in quantum fluid theory, the small momenta approximation can be applied. |22| This means neglecting 
the quadratic terms in u, v in g(u, v), which corresponds to small velocities, consistently with the SVEA. In this case 
g = (v — pu)(o;J — u! 2 )/(l + p)^" , showing that, as t — > oo, 

V pq (t,u,v) -> uj p (v)u q (u)\x pq (u + v)\ 2 2 ^ 1 / + P ^ - S(v- pu). (30) 

\U X - U3 2 \ 

Thus, in the asymptotic state the two particles travel approximately at the same velocity (ft = 1 in practical cases). 

For large propagation distances, the quantum fluid is hence a superposition of entangled pairs locked together. 
This gives a clear picture of what entanglement is to be intended for. After the generation, the two excitations are 
associated to X-shaped spatial distributions of energy that travel at the same v. Since the velocity determines the 
angular aperture (far from the origin) of the double cone which forms the X-wave.|23j a three-dimensional region 
where the two particles can be revealed simultaneously is determined. The duration of the temporal profile fixes 
the useful time window. In some sense, the paired elementary excitations are delocalized in an highly localized 3D 
X-shaped rigidly moving region. 
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CONCLUSIONS 

Using X-waves as basis, a traveling wave quantization leads to the representation of a pulsed beam as a quantum 
gas of free particles. A structured, localized, spatio-temporal distribution of energy is associated to each fundamental 
excitation. These elements furnish an alternative point of view regarding the concept of "photon" , and its localization, 
without introducing box-normalization or related ideas. 

A 3D quantum propagation is reduced to a one dimensional evolution, thus enabling simple treatments of nonlinear 
optical processes. This is mainly a consequence of the paraxial approximation, which unavoidably fixes a preferential 
direction. For this reason, this approach is well suited for describing optical experiments. It is also susceptible of 
generalizations in Lorentz covariant theories. 

During optical parametric amplification the particles composing the pulsed beam become entangled in pairs moving, 
approximately, at the same velocity. These excitations are delocalized into a 3D rigidly moving conical region, where 
they can be revealed simultaneously. The role of first order and of second order dispersion is clearly stated, while 
including diffraction. This provide a complete spatio-temporal picture of a nonlinear quantum process. 
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